integrate sqrt(1+a/x^2) dx